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Lattice QCD (LQCD) has the sign problem at real quark chemical potential. There are some 
regions with no sign problem; one is the imaginary quark chemical potential region and the others 
are the real and imaginary isospin chemical potential regions. We show that the Polyakov-loop 
extended Nambu-Jona-Lasinio (PNJL) model can reproduce LQCD data in the regions. We also 
determine the model parameters from the data and predict the QCD phase diagram in the real 
quark chemical potential region. 
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1. Introduction 



The QCD phase diagram is essential for understanding not only natural phenomena such as 
compact stars and the early universe but also laboratory experiments such as relativistic heavy- 
ion collisions. Quantitative calculations of the phase diagram from the first-principle lattice QCD 
(LQCD) have the sign problem at real quark chemical potential (/i q ). Though several approaches 
have been proposed to circumvent the difficulty, these are still far from perfection. 

As an approach complementary to the first-principle LQCD, we can consider effective mod- 
els such as the Nambu-Jona-Lasinio (NJL) model and the Polyakov-loop extended Nambu-Jona- 
Lasinio (PNJL) model. The NJL model describes the chiral symmetry breaking, but not the con- 
finement mechanism. The PNJL model is constructed so as to treat both the mechanisms. In 
the NJL-type models, the input parameters are determined at n q = 0. It is then highly nontrivial 
whether the models predict the dynamics of QCD at finite /i q properly. This should be tested from 
QCD. Fortunately, this is possible in some regions without sign problem, such as imaginary pL q , 
real and imaginary isospin chemical potential (/ii). 

In this paper, we consider two-flavor QCD and show the reliability of the PNJL model by 
comparing the model result with LQCD data in their regions. 



2. Imaginary Quark Chemical Potential 

Roberge and Weiss [jl|] found that the thermodynamic potential, £2QCD($q)> of QCD at imag- 
inary chemical potential }X q = i6 q T has a periodicity £lQCT)(Qq) = Q-QCviOq + Ink/3), showing 
that Q.QCD{6 q + 2nk/3) is transformed into Q.QCD(0 q ) by the Z3 transformation with integer k. 
This means that QCD is invariant under a combination of the Z3 transformation and a parameter 
transformation d q — >■ d q + 2nk/3. We call this combination the extended Z3 transformation. Thus, 
^QCD(6Jq) has the extended Z3 symmetry, and hence quantities invariant under the extended Z3 
transformation have the RW periodicity 

We reveal that the PNJL model has the RW periodicity [g]. The two-flavor PNJL Lagrangian [fj] 
in Euclidean spacetime is 

J&f = q{ij v D v - 74^ + m )q - G s [(qq) 2 + (qiy 5 Tq) 2 } + 0*(*[A] , 0* [A], T), (2.1) 

where q denotes the two-flavor quark field, mo does the current quark mass, and D v = d v — iA v d v o 
with the gauge field A v . In the chiral limit (mo = 0), the Lagrangian density has the exact SU (2)r x 
SU (2)l x t/(l) v x SU (3) c symmetry. The Polyakov potential t/<j> [Q] is a function of the Polyakov 
loop <& = ^tr c L with L = e' A4 ^ T and its Hermitian conjugate <&*. The PNJL thermodynamic poten- 
tial Q. in the mean field approximation (MFA) is 



Q 



2*r) 



e(p)/r+!'A6 

A=±l 



+ G s a 2 + £/<!>, (2.2) 



where e = \/p 2 + M 2 , M = mo — 2G S (7, and a = (qq). The thermodynamic potential £1 is invariant 
under the extended Z3 transformation, 

L^e- i27lk/3 L, d q ^d q + 2nk/3. (2.3) 
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Figure 1: Phase diagram in the real and imaginary jitq regions by the PNJL model with the parameter 
set [||] that reproduces the LQCD data at imaginary jitq. The points Erw and Ec are the endpoints of the 
RW transition and the first-order chiral transition respectively. The solid line denotes the first-order chiral 
transition, the dashed (dotted) line does the crossover deconfinement (chiral) transition, and the dot-dashed 
line does the RW transition. Lattice data (x) are taken from [Q], 



Therefore, Q. has the RW periodicity. 

At the present stage, the PNJL model is only a realistic effective model that possesses both 
the extended Z 3 symmetry and the chiral symmetry This property guarantees that the phase 
diagram evaluated by the PNJL model has the RW periodicity in the imaginary jU q region, and there- 
fore makes it possible to compare the PNJL result with LQCD data quantitatively in the imaginary 
/i q region. Actually, the PNJL model succeeds in reproducing the LQCD data ^ by introducing 
the vector-type four-quark interaction and the scalar-type eight-quark interaction [gj. The QCD 
phase diagram in the real /i q region is predicted by the PNJL model with the parameter set [§] that 
reproduces the LQCD data at imaginary jj, q , as shown in Fig. [j] 
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Figure 2: Phase diagrams of the deconfinement and the RW phase transition in the 9 q — T plane with 
RRW-type U$, [|]] (panel (a)) and F-type [[|] (panel (b)). The solid (dashed) line denotes the first-order 
(crossover) deconfinement transition, and the dot-dashed line does the RW transition. Point Erw is the 
endpoint of the RW transition. Lattice data (+) are taken from [||]. 

The phase diagrams of the deconfinement and the RW phase transition in the d q — T plane by 
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using the PNJL models with RRW-type [Ql and F-type [Q] are shown in Fig. ^| (a) and (b), 
respectively. Thus, the PNJL model with RRW-type reproduces LQCD data ^ at finite d q , but 
the model with F-type doesn't. In this sense, the PNJL model with RRW-type U® calculation 
is more reliable. 




Figure 3: (a) The phase structure near Ep>w with RRW-type t/<j> is magnified. The solid (dashed) line denotes 
the first-order (crossover) deconfinement transition, and the dot-dashed line does the RW transition. Points 
Erw and CP are an endpoint of the RW transition and critical endpoints, respectively, (b) T dependence of 
the chiral and Polyakov-loop susceptibilities, % a and ^<j>, at the point CP. 

The phase diagram for RRW-type near Erw is magnified in Fig. [| (a). The RW endpoint 
is first order for RRW-type U<$>, but it's second order for F-type Thus, the order of the 

deconfinement phase transition near the RW endpoint strongly depends on taken. The result of 
the PNJL calculation with RRW-type U<& is consistent with the LQCD data [^] where the order of 
the RW phase transition at E RW is first order for small quark mass. Point E RW is the triple point 
where the three first-order lines meet. Thus, there are two critical endpoints (CP) for each triple 
point ; CP is a point where the crossover and the first order lines meet. Figure [3](b) shows the chiral 
and the Polyakov loop susceptibilities, % a and as a function of T near CP. The susceptibilities 
are divergent at CP. Hence, the chiral and deconfinement transitions are second order at CP. 



3. Imaginary Isospin Chemical Potential 

LQCD has no sign problem at both real and imaginary jii. Recently, LQCD data were mea- 
sured there and also in the case where both /ii and ji q are imaginary 

In the chiral limit, QCD has the chiral SUl(2) x SUr(2) symmetry when /ii = 0. However, 
at Hi 7^ this symmetry is reduced to C/i 3 (l) x Uaj 3 (1), where t/i 3 (l) is the isospin subgroup and 
Uai 3 (1) is the axial isospin subgroup. In the case m u = m c i ^ 0, only the U\ 3 (1) symmetry survives. 
When QCD vacuum keeps the U v (l) and Uj 3 (1) symmetries, the baryon charge B = V(B) is either 
zero or integer and the isospin charge /3 = V(h) is also either zero or half-integer, where B = 
qY4q,h = qjihq and V is the volume. In the partition function Z, the baryon- and the isospin- 
charge operator appear through the form exp[V (2idih + id q B)] where ji q i = iTd q j. Therefore, Z 
has the periodicity Z(8 q , 0i) = Z(0 q , d\ +2n). In the isospin symmetric limit m u = ma, Z is invariant 
under the interchange u o d, i.e. B\ — > — d\. Hence, Z is invariant under charge conjugation, both 
d q — > — d q and di — > — d\. Furthermore we have proved that Z has the RW periodicity at B\ ^ 
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All the relations are summarized as 

Z(0 q ,0i) =Z(±0 q ,T0i) =Z(6 q , 0i + 27r) =Z(e q + 2n/3,d l ). (3.1) 

Meanwhile, if the pion condensation occurs, the Ui 3 (l) symmetry is spontaneously broken 
and the isospin charge is neither zero nor half-integer anymore. In this situation, QCD vacuum 
doesn't have the periodicities d3~l]). We have proved that the pion condensation doesn't take place 
at imaginary jii This can be understood intuitively. For real jUj, the Bose-Einstein distribution 
function has an infrared divergence at pL\ > m n /2. This induces the Bose-Einstein Condensation, 
that is, the pion condensation. For imaginary jUi, such a divergence never happens and hence no 



pion condensation occurs. As a result of this fact, Z has all the discrete symmetries (3.1). 




Figure 4: Q./T 4 , n q /T 3 and n\/T 3 as a function of 9 q and 0j. Panels (a), (b) and (c) correspond to T = 
175 MeV, while panels (d), (e) and (f) to T = 250 MeV. 

The absence of the pion condensation at imaginary pL\ is true in the PNJL model [g]. The PNJL 
thermodynamic potential at jii ^ in the MFA is 



d 3 p 



£ [3e/(p) + r £ lndet c (l+L A e - £ /(P)/ r+a ^)j + G s (a 2 + 7T 2 ) + ^,(3.2) 
/=±i L 



x=±\ 



where s± = \j (e(p) ± \h) 2 +N 2 , N = 2G s ii. The pion condensate n = (qijsTiq) is an order pa- 
rameter of the spontaneous breakings of the Ui 3 (1) symmetry. When there is no pion condensation, 
£2 is reduced to a simpler form 



d 3 p 

m~ 3 



6e(p) + T £ lndet c (l+L A e - £(p '/ T+a ^ +! '^; 
A, /=±i 



+ G s a 2 + £/ , (3.3) 



which is invariant under the extended Z3 transformation (2.3), therefore CI has the RW periodicity. 
The potential £2 has also the periodicity of d\ — > d\ + 2%. Furthermore £2 is invariant under the 
transformation, 6\ — > —6\, and also under the transformation, q — > — 8 q and L — > L T . These 
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properties guarantee that the PNJL model possesses all the symmetries in (3.1), and the model 
reproduces LQCD data [Q] qualitatively at imaginary jX\ and jj, q . 

Figure || shows Cl/T 4 , Im[w q ]/r 3 and Im[«i]/r 3 as a function of 6 q and ft in the cases of 
T = 175 and 250 MeV. Symmetries (3T) are seen in Fig. [| This result is consistent with LQCD 
[0]. If the pion condensate is nonzero, symmetries (3T) break down. Hence, the fact that 



ones 



LQCD has symmetries (3T) means that the pion condensation doesn't occur also in LQCD. As 
shown in Fig. |2| (a) for ft = 0, at temperature above Trw = 190 MeV, there is the RW phase 
transition at 6 q = tt/3 mod 2k/3, where n q = —dQ./d(iTd q ) is discontinuous. In Fig. |], T = 175 
and 250 MeV are typical temperatures below and above Trw, respectively. For any temperature, 
the RW periodicity is seen. Below Trw, these quantities are smooth at any d q and ft. In contrast, 
above Trw, O. and n\ have cusps at d q = 7i/3 mod 2n/3, while n q is discontinuous there. The 
discontinuity means the RW phase transition. Eventually, the transition occurs at d q = n/3 mod 
2k/3 when < ft < k/2 and n < ft < 2k, and at 6 q = mod 2k/3 when k/2 < ft < 3k/2 [fj. 
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Figure 5: Phase diagram of the deconfinement phase transition in the 9i — T plane. Panels (a), (b) and (c) 
correspond to q = 0, n/6 and 7C/3, respectively. The solid (dashed) line denotes the first-order (crossover) 
ttansition. The area labeled by 'RW between the two dot-dashed lines represents the region where the RW 
phase transition occurs. 

Figure |5| shows the phase diagram of the deconfinement phase transition in the ft — T plane. 
Near ft = n/2 mod n, the deconfinement phase transition is first order in all panel (a)-(c). Near 
ft = 71 mod 71, the deconfinement phase transition is first order at d q = 0, but crossover at d q = 71 / 6 
and 7t/3. The RW phase transition occurs in the area labeled by 'RW' between the two dot-dashed 
lines. 

Quantitative comparison of the PNJL model with LQCD data |7|] is made at T < T c by using the 
hadron resonance gas (HRG) model that can reproduce the LQCD data there. We have shown [^] 
that the PNJL model reproduces the LQCD data for the oscillatory patterns. For the magnitudes, 
meanwhile, the PNJL model underestimates the LQCD result. This discrepancy is understandable 
as follows. Below T c , hadronic excitations are important, but such an effect is not included in the 
MFA. By adding the hadronic correction to the PNJL model, the model agrees with the LQCD [^]. 
The HRG model works well at T < T c , but not at T > T c ; especially the HRG model doesn't 
reproduce the RW phase transition. In contrast, the PNJL model with the hadronic correction 
works both below and above T c . 



4. Real Isospin Chemical Potential 

LQCD data are available at real ji\ and jx q = [§]. The scalar-type eight-quark interaction is 
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necessary to reproduce LQCD data at imaginary ;ii q [Q]. Figure |] (a) shows the phase diagram of the 
PNJL model with the scalar-type eight-quark interaction in the }JL\ — T plane at jUq = 0. The PNJL 
model with the eight-quark interaction is also consistent with the LQCD at pi\ 7^ [10]. There is 
a tricritical point (TCP) where the first-order pion-superfluidity phase transition line is connected 
to the second-order phase transition. The critical points such as CEP and TCP are important as 
indicators of the chiral and pion-superfluidity phase transitions at compact stars and laboratory 
experiments where pL\ is nonzero generally. The TCP in the fr — T plane at jJ. q = is connected to 
the CEP in the }X q — T plane at \X\ = in the /i q — \X\ — T space JTOj], as shown in Fig. |6] (b). 



0.25 



0.2 



> 0.15 - 1 



o 



0.1 
0.05 




(a) 



0.1 



0.2 0.3 
W [GeV] 



1st order ■ 
2nd order - 
cross over o - 
cross over <1> - 
TCP 
Lattice ^ 
Lattice ^ 

0.4 



(b) 

T[GeV] 

0.2 r C 



0.1 



0.5 



A 




li q [GeV] 



m [GeV] 



Figure 6: (a) Phase diagram in the 9i — T plane at 9 q = with the eight-quark interaction. The thick- 
solid (dashed) line denotes a first-order (second-order) pion-superfluidity phase transition. The dot-dashed 
(dotted) line denotes a deconfinement (chiral) crossover transition. Lattice data are taken from [g|]. (b) Phase 
diagram in the jij — flq — T space with the eight-quark interaction. Line ABC denotes the chiral CEP, ABD 
line does the pion-superfluid TCP. The CEP and the TCP coexist on line AB. The solid (dashed) line denotes 
the first (second) order transition. 
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